Low momentum shell model effective interactions 
with all-order core polarizations 
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An all-order summation of core polarization diagrams using the low-momentum nucleon-nucleon 
interaction Vj ow _k is presented. The summation is based on the Kirson-Babu-Brown (KBB) induced 
interaction approach in which the vertex functions are obtained self consistently by solving a set of 
non-linear coupled equations. It is found that the solution of these equations is simplified by using 
Viow-k, which is energy independent, and by employing Green functions in the particle-particle and 
particle-hole channels. We have applied this approach to the sd-shell effective interactions and find 
that the results calculated to all orders using the KBB summation technique are remarkably similar 
to those of second-order perturbation theory, average differences being less than 10%. 
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Introduction. — Since the early works of Bertsch Q] and 
Kuo and Brown , the effect of core polarization in nu- 
clear physics has received much attention. Core polariza- 
tion is particularly important in the shell model effective 
interactions, where this process provides the long-range 
inter-nucleon interaction mediated by excitations of the 
core In microscopic calculations of effective inter- 
actions, core polarization has played an essential role, 
as illustrated by the familiar situation in 18 O. There 
the spectrum calculated with the bare G-matrix was too 
compressed compared with experiment, while the inclu- 
sion of core polarization had the desirable effect of both 
lowering the + ground state and raising the 4 + state, 
leading to a much improved agreement with experiment 
[HQ As pointed out by Zuker 0, the Kuo-B rown ma- 
trix elements, although developed quite some time ago, 
continue to be a highly useful shell model effective inter- 
action. It should be noted that the core polarization (CP) 
diagrams associated with the above interactions were all 
calculated to second order (in the G-matrix) in perturba- 
tion theory. But what are the effects of core polarization 
beyond second order, and how can they be calculated? 
In this Letter we would like to address these questions 
and present an all-order summation of CP diagrams for 
the sd-shell interactions. 

There have been a number of important CP studies 
beyond second order. Third-order core polarization di- 
agrams, including those with one fold, were studied in 
detail by Barrett and Kirson [|| for the sd-shell effective 
interactions. Horth- Jensen et al. [|J have carried out ex- 
tensive investigations of the third-order CP diagrams for 
the tin region. A main result of these studies is that the 
effect of the third-order diagrams is generally comparable 
to that of the second order; the former cannot be ignored 
in comparison with the latter. As is well known, high- 
order CP calculations are difficult to perform, largely be- 
cause the number of CP diagrams grows rapidly as one 
goes to higher orders in perturbation theory. The number 
of diagrams at third order is already quite large, though 



still manageable. Primarily because of this difficulty, a 
complete fourth-order calculation has never been carried 
out. It was soon realized that an order- by-order calcu- 
lation of CP diagrams beyond third order is not practi- 
cable. To fully assess the effects of core polarization to 
high order, a non-perturbative method is called for. 

In the present work, we shall use a non-perturbative 
method to carry out an all-order summation of CP dia- 
grams. Our method is based on the elegant and rigorous 
induced interaction approach of Kirson 7] and Babu and 
Brown [8j, hereafter referred to as the KBB method. In 
this formalism the vertex functions are obtained by solv- 
ing a set of self-consistent equations, thereby generating 
CP diagrams to all orders in a manner similar to the 
parquet summation |9j . Using this approach, Kirson has 
studied 18 and 18 F using a G-matrix interaction, and 
Jopko and Sprung E3 have carried out a model study 
of this approach using a separable interaction. A main 
conclusion of both studies is that when CP diagrams are 
included to all orders the effective interaction is very close 
to that give n by the bare interaction alone. In contrast, 
Sjoberg | ll| applied the Babu-Brown formalism to nu- 
clear matter and found that the inclusion of CP diagrams 
to all orders has a significant effect on the Fermi liquid 
parameters, in comparison with those given by the bare 
interaction. These conflicting results for CP studies of 
finite nuclei and infinite nuclear matter have served as a 
primary motivation for our present re-examination of the 
CP effect. 

Our application of the KBB formalism to shell-model 
effective interactions is similar to that of Kirson, but 
our treatment is different in a number of important re- 
gards. As we will discuss, the particle-core and hole- 
core coupling vertices used in the present work in- 
clude a larger class of diagrams than have been pre- 
viously studied. We will show how the inclusion of 
these diagrams is facilitated by using the recently devel- 
oped low-momentum nucleon-nucleon interaction Vjow-k 

instead of the previously used 



oped low-momentum nucJ 
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FIG. 2: Higher-order terms contributing to the vertex func- 
tions / and F, including (a) nested bubbles in bubbles and 
(b) particle-core and hole-core couplings. 



FIG. 1: Self-consistent diagramatic expansion of the ph vertex 
function, /, where £ is defined in the text. 

G-matrix. This is primarily because the G-matrix jrjlll9| 
depends on both starting energy and Pauli exclusion op- 
erator, while Mow-k depends on neither. It is noted that 
the S-wave interactions calculated from the Moszkowski- 
Scott separation method gave essentially the same matrix 
elements as Viow-k-Hjj In the subsequent Formalism sec- 
tion we will discuss these topics in greater detail, and in 
the final section we present our results together with a 
summary. Although the present calculation is restricted 
to the sci-shell effective interactions, the methods we de- 
velop can be readily expanded to other nuclear regions 
and with appropriate generalization applied to effective 
operators such as magnetic moments. 

Formalism. — The KBB induced interaction approach 
provides a very appealing way for summing up planar 
diagrams to all orders. Its fundamental requirement is 
that the irreducible vertex functions be calculated self- 
consistently. This means that any core polarization term 
contained in a vertex function must be generated self- 
consistently from the same vertex function. To illustrate 
this point, we consider the particle-hole (ph) vertex func- 
tion /. As shown in Fig. ^ / is generated by summation 
of the driving term V and core polarization terms, the 
latter being dependent on /. This then gives the self- 
consistent equation for / 



/ = V+T.g ph T,+T,g ph fg ph T,+T,g ph fg ph fg ph T,- 



(1) 



where g p h is the free ph Green function, and £ denotes 
the vertex for particle-core and hole-core coupling. The 
second-order CP diagram of Fig. 1 is the lowest-order 
term contained in "Eg^Y*. We note that for simplic- 
ity the bra and ket indices have been suppressed in 
the above equation as well as the equations that fol- 
low. For example, in Eq. the / on the LHS repre- 
sents (12 -1 |/|34 _x ), while the fifth and sixth E's on the 
RHS represent (lph^ 1 ^) and (2~ 1 \i:\p , h'- 1 4r 1 ), re- 
spectively. To further visualize the structure of /, let us 
consider £ = V in Eq. QJ: 



f = V + Vg ph V + Vg ph fg P hV + Vg ph fg P hfg P hV + 



Since / appears on both sides of this equation, it is clear 
that an iterative solution for / will yield core polarization 
diagrams to all orders, including those with "bubbles in- 
side bubbles" , like those shown in diagram (a) of Fig. [3 
For nuclear many-body calculations in general, we also 
need the particle-particle (jap) vertex function T. Like 
/, r is given by a driving term plus core polarization 
terms. Furthermore, the diagramatic representation of T 
is identical to Fig. ^ except that the hole lines 2 and 4 
are replaced by corresponding particle lines. This gives 
the self-consistent equation for T 



T = V+Y l g ph Y,+Eg ph fg ph T,+T,g ph fg ph fg ph T,- 



(3) 



To clarify our compact notation, we note that the exter- 
nal lines of the £ vertices in T are different than those 
shown in Fig. ^ The upper £ vertex, for example, now 
represents (2|£|p/i _1 4). These different £ vertices can be 
related to each other, however, via appropriate particle- 
hole transformations. 

Finally, the vertex functions / and T are coupled to- 
gether via the coupling vertex £. In the present work we 
choose 



V 



-•ph 



^ph = Vg ph V+Vg ph fg ph V+Vg ph fg ph fg ph V- 



s pp = Vg P pV + VgppFg P pV + Vg pp Tg pp rg pp V + 



(4) 



(5) 



(6) 



where g pp is the free pp Green function. 

The self-consistent vertex functions / and T are deter- 
mined from Eqs. and (PHHJl- These are similar to the 
equations used by Kirson [2j, except that our £ includes 
both ~E p h and £ pp , while the equivalent term in Kirson's 
calculations only includes £ p h, pi |21| . To see the role of 
the £ vertices, let us consider diagram (b) of Fig. Here 
the lower particle-core vertex, which contains repeated 
particle-particle interactions, belongs to £ pp , while the 
upper one, which contains repeated particle-hole inter- 
actions, belongs to Yi p h- It is, of course, necessary to 
include £ pp in order to have such CP diagrams in the 
all-order sum. Our equations are equivalent to those of 
Kirson when T, pp is set to zero, and in this case T does 



(2) not enter the calculation of /. 
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Solving the above equations for / and T may seem 
complicated, but we have found their solution can be 
simplified significantly through use of the true ph and pp 
Green functions 



G p h — 9ph + 9phfG p h 7 



Gpp 9pp 9pp^G 



vp- 



(7) 



(8) 



Using these Green functions to partially sum and regroup 
our series, the self-consistent Eqs. {Q and J3] EJ) assume 
a much simpler form 



f = V+ SG pfl S 



r = V + EG ph E, 



V + VGp h V + VGppV. 



(9) 



(10) 



(11) 



The above simplifications also aid our numerical efforts, 
and using the following iterative method we find our cou- 
pled equations can be solved rather efficiently. For the 
n th iteration, we start from and Y^ n \ to first cal- 



culate Gj^ 



and Gpp followed by S^, as seen from 
Eqs. I(7t411|). The vertex functions for the subsequent 
iteration are then obtained by taking /( n + 1 ) = V + 
-i) 



F+SW G p ™ } . The entire 



S (n) G C«) S (n) andr (»H 

iterative process begins from the initial /(°) = V+Vg ph V 
and F^ = V + Vg p hV, and typically converges after just 
a few iterations. 

In the present work, we have included folded diagrams 
to all orders. As detailed in [j^. we use this method to 
reduce the full-space nuclear many-body problem H^ n — 
E n ^ n to a model space problem H c gXm = E m x m , where 
H = H n + V, H c ff = H n + V e ff and V denotes the bare 
NN interaction. The effective interaction V g is given by 
the folded-diagram expansion 



V eS = Q-Q' Q + Q' Q Q- 



(12) 



We consider the effective interactions for valence nucleons 
in the sci-shell, and in this case Q is given by the vertex 
function T obtained earlier from the KBB equations. In 
0, the effect of higher order CP diagrams to the non- 
folded Q term was extensively studied. In the present 
work, we first calculate Q including CP diagrams to all 
orders. Then the above folded diagram series for V e g 
is summed to all orders using the Lee-Suzuki iteration 
method 0, explicitly including folded CP diagrams to 
all orders. 

For the present calculation we have chosen to use 
the low-momentum nucleon-nucleon interaction, Vi ow -k- 
Since the vertex functions / and T both depend on the 
starting energy, there would be off-energy-shell effects 
present in many CP diagrams if the G-matrix interac- 
tion were chosen. This would make the calculation very 
complicated. Viow-k, on the other hand, is energy in- 
dependent so no such difficulties are encountered. Since 



y=0.918x+0.012 



■ Matrix elements 
Regression line 



2 order matrix elements 

FIG. 3: A comparison of the second-order core polarization 
matrix elements with those of the all-order KBB calculation. 



detailed treatments ofMow-k have been given elsewhere 
El El El El El 03, here we provide only a brief de- 
scription. We define Vi ow _k through the T-matrix equiv- 
alence equations 



T(k', k, k 2 ) = V NN (k', k) + P q 2 dqV NN (k\ q) 

JO 

x T2^ T (l,k 1 k 2 ) (13) 



Tiow-k(p' ,P,P 2 ) = Vi ow ^ k (p',p) 

-P I q 2 dqViow-kip' 2 Tlow ~ k (l>P>P 2 ) 



(14) 



1 



p2 _ g2 



T{p',p,p 2 ) = Ti ow ^ k (p',p,p 2 ); (p',p)<A, (15) 

where Vnn represents some realistic NN potential 
and A is the decimation momentum beyond which the 
high- momentum components of Vnn are integrated out. 
Viow-k preserves both the deuteron binding energy and 
the low-energy scattering phase shifts of Vnn ■ Since em- 
pirical nucleon scattering phase shifts are available only 
up to the pion production threshold (Ei a i, ~ 350 MeV), 
beyond this momentum the realistic NN potentials can- 
not be uniquely determined. Accordingly, we choose 
A ps 2.0 fm -1 thereby retaining only the information 
from a given potential that is constrained by experiment. 
In fact for this A, the Mnw-k derived from various NN 
potentials [23L I24I |H [2(j are all nearly identical [lif . 
Except where noted otherwise, in our calculations we em- 
ploy the Vi ow _k derived from the CD-Bonn potential [23} . 
Results and Discussion. — As an initial study, we have 
carried out a restricted all-order CP calculation for the 
Sri-shell effective interactions. In particular, we sum only 
the TDA diagrams for the Green functions G pp and G p h, 
leaving a study of RPA diagrams to a future publication. 
Also, we use a limited shell model space consisting of 10 
orbits from Osi/2 to lpi/2 and an oscillator constant of 
7kj=14 McV. Only core excitations within this space are 
included. Vary, Sauer and Wong have pointed out 
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FIG. 4: Spectra for the 18 F (top) and 18 (bottom) systems 
calculated to different orders in perturbation theory. 



that for CP diagrams one needs to include intermediate 
states of high excitation energies (up to ~ lOTicu) in order 
for the second-order core polarization term to converge. 
In their work a G-matrix derived from the Reid soft- 
core potential was used, but our use of Vi ow _k may yield 
different results as it has greatly reduced high-momentum 
components. We plan to study this convergence problem 
for Mow-k in the near future. 

With these restrictions, we have calculated V c s from 
Eqs. (|?I- I12JI . A large number of angular momentum re- 
couplings are involved in calculating the CP diagrams. In 
this regard, we have followed closely the diagram rules in 
[28j |. In Fig. we compare the sd-shell V e fi matrix ele- 
ments calculated with all-order CP diagrams with those 
from second-order 0] CP diagrams. A least-squares fit 
was applied to the data, and it is apparent that the effect 
of including CP to all orders in our calculation is a mild 
suppression of the second-order contributions. This con- 
clusion is further born out in the calculation of the 18 O 
and 18 F spectra, the results of which are shown in Fig. 01 
Here we observe a weak suppression of the second-order 
effects in 18 O but a moderate suppression in 18 F. In the 
same figure we also observe that the spectra for different 
Viow-k derived from the CD-Bonn and Nijmegen bare 
potentials are nearly identical. 

In summary, we have presented a method based on the 
KBB induced interaction formalism for efficiently sum- 
ming core polarization diagrams to all orders in per- 
turbation theory. This summation is carried out by 
way of the KBB self-consistent equations whose solu- 
tion is significantly simplified by the use of the true pp 
and ph Green functions, and by the use of the energy- 
independent Viow-k- Although our calculation was re- 
stricted in several important aspects, we find that our 



final renormalized interaction is remarkably close to that 
of second-order perturbation theory. This is of practi- 
cal importance and a welcoming result, for it allows one 
to use the results from a second-order calculation to ap- 
proximate the contributions resulting from a large class 
of higher-order diagrams. In the future we intend both 
to expand our treatment by including additional diagra- 
matic contributions (RPA) and to use the results in cal- 
culations of other nuclear observables. 
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